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Abstract 

f"^ ' Using the 1 + 3 formulation of stationary spacetimes we show, in the 

^\ i context of gravoelectromagnetism, that the plane of the polarization of 

0^ ■ light rays passing close to a black hole undergoes a rotation. We show that 

this rotation has the same integral form as the usual Faraday effect, i.e. it 
^-l is proportional to the integral of the component of the gravomagnetic field 

I 1 along the propagation path. We apply this integral formula to calculate 

^q' the Faraday rotation induced by the Kerr and NUT spaces using the 

quasi-Maxwell form of the vacuum Einstein equations. 



PACS number(s): 04.20.-q, 04.70.Bw 
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1 Introduction 



It is a well known fact that the plane of the polarization of light rays passing 
through plasma in the presence of a magnetic field undergo a rotation which 
is called Faraday rotation (Faraday effect) [10]. One can show that a plane 
polarized wave is rotated through an angle AO given by 

A6 = -^-^ / nB {l dl (1) 

where Bu is the component of the magnetic field along the line of sight. 

It is also a well known consequence of the general relativity that light rays 
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passing a massive object are bent towards it. Several authors have considered 
the gravitational effect on the polarization of light rays by analogy with the 
Faraday effect [2,4,9]. In particular they have considered the propagation of 
electromagnetic waves in the Kerr spacetime. In [4] the authors have used the 
Walker-Penrose constant to calculate this effect for a Kerr black hole . They 
have shown that in the weak field limit the rotation angle of the plane of the 
polarization is proportional to the line-of-sight component of the black hole's 
angular momentum at the third order. In what follows we will use the Landau- 
Lifshitz 1+3 splitting of stationary spacetimes and show that the gravitational 
Faraday rotation has the same integral form as the usual Faraday effect if one 
replaces the magnetic field with the gravomagnetic field of the spacetime under 
consideration. Having found this integral form, one can use the quasi-Maxwell 
form of the vacuum Einstein equations to calculate the effect much more easily. 
In particular we show that the gravitational Faraday rotation in NUT space is 
zero, a result which needs a lot of calculation if one uses the approach based on 
the Walker-Penrose constant. 

2 1+3 formulation of stationary spacetimes (pro- 
jection formalism) 

Suppose that M. is the 4-dimensional manifold of a stationary spacetime with 
metric g b[] and p G A4 , then one can show that there is a 3-dimensional manifold 
S3 defined invariantly by the smooth map [3] 

* : M -> S 3 

where W = ^f(p) denotes the orbit of the timelike Killing vector £ t passing 
through p. The 3-space S3 is called the factor space A4/Gi, where G\ is the 
1-dimensional group of transformations generated by £ ( . Using a coordinate 
system adapted to the congruence £ t = d t we denote the projected 3-dimensional 
metric on S3 by 7 aj a (a, [3 = 1,2,3). These are the coordinates comoving with 
respect to the timelike Killing vector. One can use 7 Q( 3 to define differential 
operators on S3 in the same way that g a b defines differential operators on M.. 
For example the covariant derivative of a 3-vector A is defined as follows 

= 8 p A a + X^pA 1 

A a -f3 = d{3A a — \ 1 a pA 1 

where A"^ is the 3-dimensional Christoffcl symbol cosntructed from the compo- 
nents of j a/ 3 in the following way 

1 Note that the Roman indices run from to 3 and Greek indices from 1 to 3. 
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It has been shown that the metric of a stationary spacetime can be written in 
the following form [5] 



ds 2 = h(dx° - A a dx a ) 2 - dl 2 



(2) 



where 



and 
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h = ,g o 



dl 2 = laP dx a dx fi = {-g a p + ^^L)dx a dx 13 

9oo 

is the spatial distance written in terms of the 3-dimcnsional metric j a fj of S3. 
Using this formulation for a stationary spacetime one can write the vacuum 
Einstein equations in the following quasi-Maxwell form [6] 



div B n 







Curl E 3 = 



div V g = -(±(VhB g ) 2 +E 2 ) 
Curl (VhB g ) = 2E 3 x (VhB g ) 
pa p = E a;i3 + ([y/hB2)(y/hB$) - (VhB g ) 2 ^ 
where the gravoelectromagnctic fields arc 
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(3) 
(4) 
(5a) 

(56) 
(6) 

(7) 
(8) 

and P af3 is the 3-dimensional Ricci tensor constructed from the metric -f a p. It is 
attractive to regard the combination \fhB g , appearing in the above equations, as 
the gravitational analogue of the magnetic intensity field H and denoting it with 
H g . In this way one may think of the last term in equation (5b) as an energy 
current corresponding to a Poynting vector flux of gravitational field energy. 
Note that all operations in these equations are defined in the 3-dimcnsional 
space with metric y a /3- Using the timelike Killing vector of the spacetime one 
can define the above gravoelectromagnetic fields in the following covariant forms 



B g = Curl A. 



E, 



1 (ego) 

2 Id 2 



1*1 = h 1 ' 2 



where e b ^ d is the 4-dimensional antisymmetric tensor and "; " denotes covariant 
differentiation. 
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3 Derivation of the gravitational Faraday rota- 
tion 



We use the analogy with the flat spacetime and take the plane of the polarization 
of an electromagnetic wave to consist of two 3-vectors k and f , the wave vector 
and the polarization vector respectively. The 4-vectors corresponding to these 
two 3-vectors have the following relations 

k a k a = k a f a = /°/a = l a = 0,1,2,3 (9) 

Both of these 4-vectors are parallely transported along null geodesies [7] i.e. 

dk a 

v k k a = — + r a mn k n k m = o 

V k f a = ^+T a mn f n k m = 

where A is an affine parameter varying along the ray. Employing an orthogonal 
decomposition based on the adapted coordinates, the above 3-vectors defined on 
the 3-space S3 can be taken to be equivalent to the contravariant components of 
k a and f a i.e. k = ^k a = ^k a and f = ( 3 )/ a = ^ f a [8]. One should note 
that the covariant counterparts of these 3-vectors are not the spatial components 
of the covariant 4-vectors k a and f a but 

and 

(3) //3 = 7 Q/3 (3) r = {i) fp + kg P - 

From equation (9) one can see that the polarization vector is known up to a 
constant multiple of the wave vector i.e. both f a and f' a = f a + Ck a satisfy 
equation (9). This shows that there is a kind of gauge freedom in choosing 
/ which enables one to put fo = without loss of generality and in which 
case ( 3 ) fp = < - 4 ' ) //3-Q Applying the above decomposition the evolution of the 
3-vectors k and f along the ray is given by the spatial components of the parallel 
transport equations i.e. 

f)k a 

— + T« in k n k m = (10) 

^ + r«„rfc m = o (11) 

Now we try to write these two equations in terms of the 3-dimensional quantities 
defined on £3. From equation (10) we have 

r)h a 



2 This choice corresponds to C = — and makes / orthogonal to the time lines. 
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to calculate this we need the following components of the christoffel symbol [5] 

V a — -h- a 
1 oo — ~ n > 
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^ = ~{9°b-9f)-\90hf (12) 

r^ 7 = a^ 7 + |^(fl- Q - <£) + 57 « - <#)] + \m^? 

Substituting the above equations into equation (10) we have 

3 V k k=-fc (B 9 xk) + ^E 9 (13) 



where 



is the 3-dimensional analogue of the parallel transport of k along itself in £3, 
and we have used equations (7) and (8) and the fact that 

fc° = ^ + g.k (14) 

Now in the same way we write equation (11) in terms of 3-dimensional quantities 
defined on S 3 

a fa 

Substituting from equations (12) and using the gauge in which /o = we have 

3 V k f = ~k {B g xf) (15) 

To be able to interpret equations (13) and (15) we need to perform another 
calculation. Using the facts that f .f = 1 and f .k = one can show that F] 



-fc (B 9 x k) - -° (E ff .f)f - -fc (B s .f)(f x k) (16) 



1 

2" uv -»"-' h"~ y "'' 2 

Using the above equation one can write equations (13) and (15) in the following 
forms 

3 V k k = Lxk+(E s .k)k (17) 



3 We have also used the relation 

k 2 



ifc (B 9 xf).k = --&E 8 .f 
2 h 



which can be found by considering the fact that 

3 V k (k.f) = ( 3 V k k).f + k.( 3 V k f) = 0. 
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3 V k f = Lxf (18) 

where 

L = -h Q [B g - i(B s .f)f + |k|(E g .n)f] (19) 

where n = f x k is a unit vector in the polarization planeQ. If we had only 
the second term in the RHS of equation (17) that would have meant, by com- 
parison with the 4-dimensional definition of the parallel transport, that the 
3-dimensional vecor k is parallely transported along the projection of the null 
geodesic in E 3 space. But the appearance of the first term shows that k has 
also been rotated by an angular velocity L. The same rotation happens to the 
polarization vector f as can be seen from equation (18). Therefore the combi- 
nation of these two equations leads us to the fact that the polarization plane 
has rotated by angular velocity L along the projected null geodesic. What we 
are interested in is the angle of rotation around the tangent vector k along the 
path between the source and the observer. So we have 

robs. 

n= L.kdA 

J sou. 

Substituting for L from (19) we have 

1 pobs. 

fl = — / k B g .k dX (19a) 

2 J sou. 

where we used the fact that f .k = 0, which follows from equation (9) and the 
gauge freedom which allows the choice /o = 0. Now combining the following 
two equations 

k = g 0a k a = h(k° - g a k a ) 
k a k a = = h(k° - g a k a f - jafj^k 13 = 

we have 

la pk a k? = 



k 2 



or equivalently in terms of k a = — 



finally upon substitution of (19b) in (19a) and putting h.dl = dl we find 

i pobs. ^ 

n = — / VhB g .dl (20) 

^ J sou. 

which has the same integral form as equation (1) i.e. the gravitational Faraday 
rotation is proportional to the integral of the component of the gravomagnetic 

4 The same relations were also found by Fayos and Llosa [2] where they have not used the 
unit vector n that we have used here. 
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field along the propagation path. But their main difference is the fact that 
the gravitational Faraday rotation, given by (20), is a purely geometrical effect 
while the usual Faraday effect, equation (1), depends on the frequency of the 
light ray. In the next two sections we will apply this formula to the cases of 
NUT and Kerr black holes. 



4 Gravitational Faraday rotation in NUT space 

There is no gravitational Faraday rotation induced by NUT space and the reason 
is as follows. Take a closed path C around the NUT hole which consists of two 
paths (see figure 1). Path 1, a null geodesic which passes close to the black hole 
and path 2 so far away that the effect of the gravitational field (including the 
Faraday rotation) on the light rays is negligible (another reason that on path 2 
there is no gravitational Faraday rotation is the fact that B 9 — > as r — > oo.). 
Now using the Stokes theorem, one can write equation (20) in the following form 

and using equation (5b) we have 

-\ J{VhB g ).d\ - \ J{VhB g ).d\ = - J(V g x VhB g ).dS (21) 

The second term in the LHS of the above equation is zero by costruction. 
On the other hand for the NUT space we have^j 

E s = -ld r (lnf(r))r 



n = --6 (VhB g ).dl = -- / V x (V/iB 9 ).dS 



and 



_ 2Z/(r)s 
a — „2 



Which together show that the RHS of equation (21) is also zero and therefore 

/ ^B s ).dl = (22) 



1 

~2 

i.e. there is no Faraday effect on the light rays passing a NUT black hole close 
by. One can show that the same result can be obtained for the NUT space 
using the approach based on the Walker-Penrose constant. In this case one 
needs to take into account the simplifying fact that all the geodesies in NUT 
space including the null ones lie on spatial cones [6] . 

5 We have used the following form of the NUT metric 

ds 2 = f(r)(dt - 2lcos8d4) 2 - /(r)" 1 *- 2 - (r 2 + l 2 )(d9 2 + sin 2 6d<p 2 ), 

where 

J( ' (r 2 + l 2 ) 
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5 Gravitational Faraday rotation in Kerr metric 



Faraday effect in Kerr metric has already been studied and it has been shown 
that despite previous claims [9], when light ray passes through the vacuum 
region outside rotating matter its polarization plane rotates [4]. In this section 
we will consider Two different cases: 

1- When the orbit lies in the equatorial plane i.e. for 9 = 0. 

2- A more general orbit which intersects the equatorial plane and is symmetric 
about it. 



5.1 Orbits in the equatorial plane 

In this case using the definitions of E g and B ff and equation (21) one can see 
that the gravitational analogue of Poynting vector defined by E 5 x VhB g has 
only one component along the <j> direction and therefore is normal to the plane 
of the orbit, which in turn leads to the fact that in this special case there is no 
gravitaional faraday effect on light rays. 



5.2 A symmetric orbit about the equatorial plane 



In this case we need to find the orbit and we will see that one just needs to find 
the orbit in the zcroth order in a/r and m/r (i.e. straight line approximation) 
which is done in appendix A. 

Writing the Kerr metric in form (2) in Boyer-Lindquist coordinates one can see 
that 

2aTOrsin 2 # 

A = A d 



2mr — p 2 



from which we have 



and 



B„ 



2amrsm29[2mr 



^y{2mr — r 2 — a 2 cos9 ) 

2a77isin 2 #(a 2 cos# 2 — r 2 ) 
y/~f(2mr — r 2 — a 2 cos 2 9) 2 



where 



7 = det^ap and p 2 = r 2 + a 2 cos 2 1 



Using the definition of the gravoelectric field given in (7) we have 



and 



Am(a 2 cos 2 # - r 2 
p 4 (p 2 -2mr) 

^ g rma 2 sm29 
9 ~ p A {p 2 - 2mr) 
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where A = r 2 + a 2 — 2mr. Substituting the above fields in equation (21) and 
putting \i — cos9 we have 



O = - f(E g x VhBg)^^ = 2am 2 ^ j 

J s J —f-i J r 



^ - 2 (23) 

.(n) (r 2 + a 2 y? — 2mr) 

where r or b.(ii) is the equation (of the projection ) of the orbit in the (r, 0) plane. 
To find the lowest order Faraday effect we calculate the above integral neglecting 
the a 2 /r 2 and m/r terms in which case we have 

/•Mo poo ^ ^ pHo ^ 

£1 = 2am 2 / / — drdyU = — -am 2 / ^-d/x 

JO forb 



Jr orh . 0*) r4 3 Jo r or b 3 



Using the (r, 6*) equation of the orbit given in appendix A one can calculate 
the above integral which gives 



Q = -\am 2 f (1 - !f^ M 2 )3/2 rf/i = (i/4)^ cos0o ^l 



This expression is of the third order r ° m which is of the same order as the 
result given in [4]. 

Discussion 

We have shown that using the 1+3 formulation of stationary spacetimes one can 
cast the gravitational Faraday rotation in exactly the same mathematical form 
as the usual Faraday effect i.e the gravitational Faraday rotation is proportional 
to the gravomagneic field of the spacetime along the propagation path of the 
light ray. One should note that the origins of these two effects are completely 
different. The usual Faraday effect originates from the interactions between the 
electrons in plasma in one hand with the electromagnetic field of the light ray 
and the external magnetic field on the other hand and therefore depends on the 
frequency of the light ray. But the gravitational Faraday effect is a purely geo- 
metrical one originating from the structure of the spacetime under consideration 
and it is normally attributed to the behaviour of the reference frames outside 
the spacetime of a rotating body. Using the quasi-Maxwell form of the vacuum 
Einstein equations we showed that there is an easy way to calculate the effect 
by transforming the line integral of B s to a surface integral of the gravitational 
analogue of the Poynting vector. More importantly the order of the effect can 
be seen without going through the detailed calculation (as in eqation (23) for 
the Kerr case) and in some cases like NUT space just a simple observation re- 
veals that there is no effect at all. For gravitational waves of small amplitude 
propagating in a curved background, one can develope the geometric optics in 
such a way that the wave and polarization 4-vectors satisfy the same relations 
as given by equations (2) [7]. So it can easily be seen that all the main relations 
that we have found for the gravitational Faraday rotation of light rays are also 
applicable to gravitational waves of small amplitude. 
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Appendix A 

The equation governing the projection of the orbit in the (r, 9) plane for Kerr 
metric is given by [1] 

dr 



^Jr 4 + (a 2 — £ 2 — rf)r 2 + 2m[r] + (£ — a) 2 ]r — a 2 rj 

de 



\fr) + a 2 cos 2 9 — £ 2 cot 2 9 



(Al) 



where £ and r\ are the constants of the motion and we choose the case in which 
T) > 0, which corresponds to the null geodesies which intersect the equatorial 
plane and are symmetric about it [1]. We perform the above integrations for 
the case when a/r <C 1 and m/r <C 1 i.e. for weak deflections and indeed as we 
will see for a case in which there is no deflection in the (r, 9) plane. First we 
evaluate the LHS of the above equation which can be written in the following 
form (after discarding the small terms) 

/*/-> v/T^-7?=< !/„„„,>„»,,,„„/.•) (.42) 



where r m i n = ^/£ 2 + r\ is the leading term (in the expansion) of the largest root 
of r 4 + (a 2 — £ 2 — rj)r 2 + 2m[rj + (£ — a) 2 ]r — a 2r q = for small deflection [4]. 
Now we evaluate the RHS of the equation (Al) in the same limit. This integral 
can be written in the following form 



RHS = 



d/i 



^Jrj + n 2 (a 2 -£ 2 -rj) -a 2 ^ 4 



Now using the fact that a/r min <C 1 and r min = ^/£ 2 + 77 we can approxiamte 
and evaluate the above integral as follows 

RHS = - [ d ^ = = -(l/r min )arcsin(^V r 2 min /v) ( A 3) 

J VV- V r min 

Equating the equations (A2) and (A3) we will have 

r ° rb = h ( l mm , (j44) 
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which is the projection of the orbit in the (r, 9) plane for small deflections and in 
this case in fact no deflection because there is no term depending on m or a. As 
one can see r — > oo when cos0 = ±^^- where plus and minus signs correspond 
to the position angles 8 Q and 9 S of the observer and the source respectively. 
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Observer 



Figure 1: The NUT hole and a closed path C around it. 



13 



